Upper Bound on the Dark Matter Total Annihilation Cross Section 
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We consider dark matter annihilation into Standard Model particles and show that the least de- 
tectable final states, namely neutrinos, define an upper bound on the total cross section. Calculating 
the cosmic diffuse neutrino signal, and comparing it to the measured terrestrial atmospheric neutrino 
background, we derive a strong and general bound. This can be evaded if the annihilation products 
are dominantly new and truly invisible particles. Our bound is much stronger than the unitarity 
bound at the most interesting masses, shows that dark matter halos cannot be significantly modified 
by annihilations, and can be improved by a factor of 10-100 with existing neutrino experiments. 

PACS numbers: 95.35.+d, 98.62.Gq, 98.70.Vc, 95.85.Ry 



The self-annihilation cross section is a fundamental 
property of dark matter. For thermal relics, it sets the 
dark matter mass density, Qdm ~ 0.3, and in these and 
more general non-thermal scenarios, also the annihila- 
tion rate in gravitationally-collapsed dark matter halos 
today Q. How large can the dark matter annihilation 
cross section be? There are two general constraints that 
bound the rate of dark matter disappearance. (Through- 
out, we mean the cross section averaged over the halo 



velocity distribution, i.e., ((Tav), where v r 
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The first is the unitarity bound, developed for the early 
universe case by Griest and Kamionkowski and for 
the late-universe halo case by Hui [3j. In the plane of 
(o~av) and dark matter mass m x , this allows only the 
region below a line (ctav) ~ V m x (this will be made more 
precise below). The second is provided by the model 
of Kaplinghat, Knox, and Turner (KKT) 4], in which 
significant dark matter annihilation is invoked to resolve 
a conflict between predicted (sharp cusps) and observed 
(flat cores) halo profiles. Since this tension may have 
been relaxed [l[, we reinterpret this type of model as 
an upper bound, allowing only the region below a line 
(<jav) ~ rn x . That the KKT model requires (cta v ) values 
> 10 7 times larger than the natural scale for a thermal 
relic highlights the weakness of the unitarity bound in 
the interesting GeV range. However, there have been no 
other strong and general bounds to improve upon these. 

While these bound the disappearance rate of dark mat- 
ter, they say nothing about the appearance rate of annihi- 
lation products, instead assuming that they can be made 
undetectable. To evade astrophysical limits, the branch- 
ing ratios to specific final states can be adjusted in model- 
dependent ways. However, a model-independent fact is 
that the branching ratios for all final states must sum to 
100%. A reasonable assumption is that these final states 
are Standard Model (SM) particles. We show that the 
most difficult SM final state to detect is neutrinos; but 
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FIG. 1: Annihilation of dark matter into SM final states. 
Since all final states except neutrinos produce gamma rays 
(see text), we can bound the total cross section from the neu- 
trino signal limit, i.e., assuming Br ("Invisible") ~ 100%. 



that surprisingly strong flux limits can be simply derived 
from recent high-statistics data; and that we may inter- 
pret these as bounding all SM final states, and hence the 
dark matter total annihilation cross section. See Fig. [TJ 

If dark matter is not its own antiparticle and if there is 
a large particle-antiparticle asymmetry, then annihilation 
could be prohibited, making all bounds inapplicable or 
irrelevant. Our bound can be evaded if the final states 
are dominantly new and truly invisible non-SM particles, 
in which case all dark matter annihilation searches will 
be more challenging; we quantify an upper bound on the 
branching ratio to SM final states below. 

Probing Dark Matter Disappearance. — For dark 
matter that is a thermal relic, the cross section required 
to ensure Q, DM ~ 0.3 is (<r A v) - 3 x 10" 26 cm 3 s" 1 0. 
KKT discussed several models in which the dark matter 
is not a thermal relic, e.g., it might have acquired mass 
only in the late universe, or have been produced through 
the late decays of heavier particles As emphasized 
in Refs. 0, 0] , it is interesting to directly ask how large 
the annihilation cross section could be in halos today, 
irrespective of possible early- universe constraints. 

Unitarity sets a general upper bound on (ctavY, and 
can only be evaded in certain unusual cases @, S @] • In 
the low-velocity limit where the cross section is assumed 
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to be s-wave dominated, (o~av) < 47r/m 2 u, or 



(a A v) < 1.5 x 10 
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cm" 



GeV 



300 km/s 



(1) 



In the KKT model, the required cross section to suffi- 
ciently distort the dark matter profiles of galaxies is 

,3 



(CAI^KKT ~ 3 X 10 



-19 



GeV 



(2) 



(Similar effects are attained via elastic interactions [6[ . A 
large self-annihilation cross section implies a large elas- 
tic self-scattering cross section, but not vice- versa 
Hui Q showed that unitarity restricts the KKT model 
to relatively small masses; for v rms — 300 km/s, m x < 
80 GeV. There have been no other model-independent 
methods to constrain the KKT model. We argue next 
that dark matter disappearance must be accompanied 
by the appearance of something, and the bound on the 
weakest final state bounds all of them. The appearance 
rate bounds {<Ja v ) directly, independent of which partial 
waves dominate a a, i-e., its u-dependence. 

Revealing Neutrino Appearance. — We assume 
that annihilation proceeds to SM particles, and express 
the cross section in terms of branching ratios to "visi- 
ble" and "invisible" final states, such as gamma rays and 
neutrinos, respectively, as in Fig. [TJ If the branching ra- 
tio to a specific final state were known, then a bound on 
that appearance rate would yield a bound on the total 
cross section, inversely proportional to this branching ra- 
tio. However, the branching ratios are model-dependent, 
and any specific one can be made very small, making 
that bound on (a^w) very weak, e.g., for m x = 1 GeV, 
KKT require Br(j) < 10~ 10 to allow their total cross 
section 4]. Note that gamma-ray data constrain only 
the product (<jav) Br(y), and the bounds vary with m x . 

KKT 0] and Hui Q assume invisible but unspecified 
final states. It is clear that most SM final states pro- 
duce gamma rays. Quarks and gluons hadronize, pro- 
ducing pions, where tt° — > 77; the decays of weak bosons 
and tau leptons also produce 7r°. The stable final state 
e + e~ is not invisible, since it produces gamma rays ei- 
ther through electromagnetic radiative corrections [3] or 
energy loss processes |8|; the final state [i + [i~ produces 
e + e~ by its decays. Thus the only possible "invisible" 
SM final states are neutrinos. 

Of final states with neutrinos, we focus on vv. Sim- 
ilar bounds could be derived for vvvv, but we as- 
sume that these are suppressed and/or that the Rube 
Goldberg-ish Feynman diagrams required would con- 
tain charged particles, and hence gamma rays through 
(model-dependent) radiative corrections. Due to elec- 
troweak bremsstrahlung, final-state neutrinos are in- 
evitably accompanied by weak bosons and hence gamma 
rays, primarily with E 1 ~ m T /2; however, these gamma- 
ray constraints on (<jav ) are weaker than or comparable 
to what we obtain directly with neutrinos [9( . 



To derive our bound on the total annihilation cross 
section, we assume Br(vv) ~ 100%. This is not an as- 
sumption about realistic outcomes, but it is the right way 
to derive the most conservative upper bound for SM final 
states. Why is this a bound on the total cross section, and 
not just on the partial cross section to neutrinos? Sup- 
pose that Br(vv) were reduced enough that the 1 / Br(vv) 
correction for an impure final state was necessary; at our 
factor-two precision, this occurs when another SM final 
state has a comparable branching ratio. For the total 
cross section set by the neutrino bound, any other pure 
final state would be more strongly constrained, thus mak- 
ing this cross section disallowed for all final states in the 
SM. Therefore, while setting this bound using neutrinos 
can be too conservative, it can never overreach. 

Cosmic Diffuse Neutrinos: Signal. — The most 
direct approach to bound the XX ^ vv cross section is 
to use the cosmic diffuse neutrino flux from dark matter 
annihilations in all halos in the universe as the signal. 
Since this is isotropic and time-independent, it is chal- 
lenging to detect above the background caused by the 
atmospheric neutrino flux. A complementary approach 
uses the Milky Way signals, which have somewhat differ- 
ent uncertainties on the predictions and data [Io| . The 
data to test the diffuse signal are available in the full 
energy range now, but this is not yet true for all the di- 
rectional signals. While the latter will likely be stronger 
eventually, going beyond our rough estimates will require 
proper experimental analyses. 

The cosmic diffuse signal from XX ~> vv annihilations 
depends on the radial density profile of each dark matter 
halo, the halo mass function (the relative weighting of ha- 
los of different masses) , and how those halos evolve with 
redshift. We follow the calculations of Ullio et 
see also [l3j]. The signal spectrum is 
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dz 



A 2 {z)dN iy {E') 
h(z) 



dE' 



dE 2 AttHq ml 

(3) 

where the 1/2 is for assuming x is its own antiparticlc, 
Hq = 100 ftkms -1 Mpc -1 is the Hubble parameter, and 
Qdm is the dark matter density in units of the critical 
density p C rit- We assume a flat universe, with flnM = 0.3, 
fl A = 0.7, h = 0.7, and h(z) = [(1 + zfVl DM + ft A ] 1/2 - 

The factor A 2 (z) accounts for the increase in density 
due to the clustering of dark matter in halos, defined so 
that A 2 = 1 corresponds to all dark matter being at its 
average density in the universe today. The concentration 
of halos and thus the value of A 2 evolves with redshift. 
(Note that we have absorbed a factor of (1 + z) 3 into the 
definition of A 2 , as in Ref. 12j.) However, to collect most 
of the signal, we only need neutrino energies near m x , and 
hence will be sensitive only to modest redshifts where it is 
accurate to take A 2 (z) ~ A 2 (0) [11]. Note that only A 2 
matters, and not its individual factors. The value of A 2 
does depend on the halo profile chosen. We adopt A 2 = 



3 



-a 



O -10 
W -12 



M 
O 



10 



i r 



j L 



data 



assumed 



-1 1 2 3 4 5 6 
logE [GeV] 

FIG. 2: Upper: Diffuse vv annihilation signal for m x — 
10, 10 3 , and 10 s GeV, added to the atmospheric background, 
both as (y^ + Vfj) and versus neutrino energy. As noted, the 
signals are most accurate for E v > m x /3. Lower: Ratio of 
this sum and background. The (oav) values at each example 
m x are chosen to be detectable by our conservative criteria; 
the data and assumed uncertainty scales are also indicated. 



2 x 10 5 for Navarro, Frenk, and White (NFW) halos with 
moderate assumptions about the halo mass distribution. 
For cuspier Moore profiles, A 2 could be ~ 10 times larger, 
while for flatter Kravtsov profiles, it could be ~ 2 times 
smaller; see the discussions in Ref. (lo| . 

For XX vv, the source spectrum dN u /dE' is a delta 
function; neutrinos produced with energy E' are red- 
shifted to the observed energy E = E' /(l + z), i.e., 



dN„(E') 
dE' 



= -6(m v -E') = — 5 
3 v x ' 3E 



E 



- 1 



(4) 



where we have accounted for 2 neutrinos per annihilation, 
equally divided among 3 flavors. (Note that has a 
large fraction in every neutrino mass eigenstate, so any 
initial mix of mass or flavor eigenstates would be close 
to this.) In Fig. [21 we show example dark matter signals 
compared to the atmospheric neutrino background. 

Using the neutrino signal, we can also derive a con- 
straint from the relativistic energy density. Requir- 
ing f2 ra d < 0.2 at low redshift [l4j leads to (<jav) < 
10~ 17 x (m x /GeV) cm 3 /s. While this bound applies to 
any light final state including non-SM particles such as 
purely sterile neutrinos, it is weak, and would require 
even greater halo modifications than the KKT model. 

Cosmic Diffuse Neutrinos: Backgrounds. — How 
large of a neutrino signal is allowed by present data? As 
shown in Fig. [21 the signal spectrum is sharply peaked. 
To be insensitive to the spectrum shape, i.e., the red- 
shift evolution, we define the signal as integrated over 



a bin of width Alog^E'j/ = 0.5, just below E v = m x , 
(i.e., we consider z < 2, so we can ignore the tails at 
low energy.) To be detectable, we require that the sig- 
nal be 100% as large as the angle- averaged atmospheric 
neutrino (y^ + v^) background, integrated in the same 
way. Signal and background are both somewhat smeared 
from received neutrino energy to detected energy, but 
well within this bin. This conservative approach allows 
us to simply derive the flux and annihilation cross section 
constraints over the very wide mass range 0.1-10 5 GeV. 
Our model predicts equal fluxes of (y e + T> e ), (y^ + v^), 
and (z/ r + P T ), any of which can be used to derive bounds. 

The atmospheric neutrino (y^ + u^) spectra as a func- 
tion of neutrino energy have been derived from data 
from the the Frejus (0.25-10 4 GeV, in 9 binsljp and 
AMANDA (1.3xl0 3 -3.0xl0 5 GeV, in 10 bins) jnfdetec- 
tors. Neutrino attenuation in Earth will only be signifi- 
cant above 10 5 GeV. The agreement with theoretical pre- 
dictions against upward fluctuations in the data is very 
good, well below the 100% uncertainty that we adopted. 
These spectra were derived from neutrino-induced muon 
data by a regularized unfolding technique, which might 
miss a narrow signal. 

We thus considered the data in more detail, find- 
ing that for E v — 0.1-10 4 GeV, such a signal is defi- 
nitely excluded, especially using both the (y^ + v^) and 
(v e + D e ) signals. The most useful data are from the 
Super-Kamiokande detector. In Ref. [I?]], visible-energy 
spectra for each of e-like and ^t-like events from 0.1-100 
GeV are given in 4 log-spaced bins per decade. The 
agreement with predictions including neutrino oscilla- 
tions is excellent; the moderate exce ptio ns in some of the 
highest-energy bins are explainable 17] . Neutrinos with 
E v ~ 10-10 3 GeV are probed by the count rates (no spec- 
tra) of upward throughgoing muons 17| , and similarl y fo r 
E v ~ 10 2 -10 4 GeV and upward showering muons [18j |: 
both are also in excellent agreement with predictions. 

Dedicated analyses of the measured data could improve 
the signal sensitivity by a factor 10-100, depending on 
the energy range. First, using the sharp feature in the 
spectrum at m x ; Fig. [2] shows that while the signal is 
comparable to the background when integrated over an 
energy bin, in the endpoint region it is much larger. Sec- 
ond, the uncertainties below 10 GeV are actually below 
10%, and apply to narrower bins in energy than we as- 
sumed [l3|. Third, by 10 GeV, the (i/ e + D e ) to ( v„, + v^) 
background flux ratio is 1/3 and rapidly falling [19|; in 
addition, the (y e + u e ) flux is strongly peaked at the hori- 
zon (l9j . while the dark matter signal is isotropic. Fourth, 
detailed analyses of Super-Kamiokande and AMANDA 
upward throughgoing and upward showering muon data 
should be more sensitive than the simple count rates we 
used. 

Conclusions. — We have shown that the dark mat- 
ter total annihilation cross section in the late universe, 
i.e., the dark matter disappearance rate, can be directly 
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FIG. 3: Upper bounds on the dark matter total annihilation 
cross section in galaxy halos as a function of the dark matter 
mass, calculated as discussed in the text. 



and generally bounded by the least detectable SM states, 
i.e., the neutrino appearance rate. This can be simply 
and robustly constrained by comparing the diffuse signal 
from all dark matter halos to the terrestrial atmospheric 
neutrino background. Our bound on {uav) is shown in 
Fig. El Over a large range in m x , it is much stronger than 
the unitarity bound of Hui Q . ft strongly rules out the 
proposal of Kaplinghat, Knox, and Turner Q to mod- 
ify dark matter halos by annihilation. Our bound can 
be evaded with truly invisible non-SM final states. For 
a cross section above our bound, its ratio to our bound 
yields an upper limit on the branching ratio to SM final 
states required to invoke that large of a cross section. 

Annihilation flattens halo cusps to cores of density 
Pa ~ m x/(( (7 Av)HQ 1 ) Our bound implies that for 
all m x > 0.1 GeV, this density is pa > 5 x 10 3 GeV/cm 3 , 
which only occurs at radii < 1 pc in the Milky Way for 
an NFW profile. Annihilation should thus have minimal 
effects on galactic halos. 

Detailed analyses by the Super-Kamiokande and 
AMANDA Collaborations should be able to improve our 
bound by a factor 10-100 over the whole mass range. 
Halo substructure or mini-spikes around intermediate- 
mass black holes could increase the signal by orders of 
magnitude [13, [20j . The sensitivity could thus become 
close to the natural scale for thermal relics, making it a 
new tool for testing even standard scenarios. 
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